Quantum jumps and entropy production 
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The irreversible motion of an open quantum system can be represented through an ensemble of 
state vectors following a stochastic dynamics with piecewise deterministic paths. It is shown that 
this representation leads to a natural definition of the rate of quantum entropy production. The 
entropy production rate is expressed in terms of the von Neumann entropy and of the numbers of 
quantum jumps corresponding to the various decay channels of the open system. The proof of the 
positivity and of the convexity of the entropy production rate is given. Monte Carlo simulations 
of the stochastic dynamics of a driven qubit and of a A-configuration involving a dark state are 
performed in order to illustrate the general theory. 

PACS numbers: 03.65.Yz, 02.70.Ss, 42.50.Lc 

Keywords: open quantum systems, quantum entropy, quantum jumps, stochastic wave function method 



I. INTRODUCTION 

Dynamical processes in dissipative open systems in- 
volve, in general, irreversible transformations of non- 
equilibrium states Q, @] . An important physical quantity 
of the theory of non-equilibrium processes is the entropy 
production, that is the rate at which entropy is produced 
as a result of irreversible processes Q. In the classi- 
cal theory the entropy production provides an important 
Lyapunov functional of non-equilibrium stationary states 
and leads to the principle of minimal entropy production. 
It further enables the determination of non-equilibrium 
transport properties and Onsager coefficients (for a re- 
cent discussion of the classical theory see, e. g. 0,0) 

In the quantum theory of open systems |(| the entropy 
production is usually related to the negative derivative 
of the relative entropy with respect to an invariant state. 
On the ground of the weak coupling master equation for 
the open system's density matrix the corresponding ex- 
pressions can again be shown to have useful convexity 
properties and to satisfy a minimal entropy production 
principle 0, as well as a kind of quantum Onsager rela- 
tions The concept of relative entropy is also useful 
in quantum information theor y |9| , e. g. in the quan- 
tification of entangleme nt JlOl [llj and of the success of 
quantum teleportation |12j . A general account of the 
theory of quantum entropy may be found in [l3|. 

During the past decade there has been considerable 
interest in the application of stochastic wave function 
methods to the analysis and simulation of dissipative 
open quantum systems. By contrast to the conven- 
tional description in terms of the reduced density ma- 
trix, in these methods one considers an ensemble of pure 
states whose time-evolution is governed by a stochas- 
tic process in the underlying Hilbert space Ti of the 
open system. In a particularly interesting method of 
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this kind, developed by several authors 0,0,0, the 
wave function represents a piecewise deterministic pro- 
cess in Hilbert space. This means that the stochastic 
evolution consists of smooth, deterministic parts which 
are interrupted by instantaneous quantum jumps. Re- 
cently, it has been shown that the stochastic wave func- 
tion method can be generalized to enable the treatment 
of non-Markovian quantum processes [l7j and of Bosonic 
jl8j and of Fermionic |19| many-body systems by means 
of an appropriate stochastic representation in the dou- 
bled Hilbert space H®H. 

The aim of the present paper is to connect the stochas- 
tic state vector dynamics to the production of quantum 
entropy. To this end, we apply the entropy balance equa- 
tion of non-equilibrium thermodynamics to the von Neu- 
mann entropy of the open system. A quantum jump 
of the stochastic dynamics is linked to the exchange of 
an entropy quantum between system and environment. 
This enables one to express the entropy flux in terms of 
the random processes which count the numbers of quan- 
tum jumps belonging to the various decay modes of the 
open system. In this way, a stochastic representation of 
the production rate of quantum entropy is derived from 
the piecewise deterministic state vector evolution. It will 
be demonstrated that the obtained expression for the 
entropy production provides a non-negative and convex 
functional on the state space of the open system. Fur- 
thermore, the stochastic representation is shown to lead 
to a Monte Carlo simulation method for the determina- 
tion of the quantum entropy production in open systems. 

The paper is organized as follows. In Sec. II we 
briefly recall the dynamical description of open quantum 
systems in terms of quantum master equations and of 
stochastic differential equations. Sec. Ill contains the 
construction of the stochastic representation of the en- 
tropy production rate as well as the proof of its posi- 
tivity and convexity. Finally, a discussion of the results 
and examples of Monte Carlo simulations of the entropy 
production are given in Sec. IV. 



2 



II. QUANTUM MASTER EQUATIONS AND 
PIECEWISE DETERMINISTIC PROCESSES 

We consider an open quantum system with reduced 
density matrix p = p(t). The open system is coupled to 
a heat bath, that is to an external reservoir in a thermal 
equilibrium state at temperature T. The free evolution of 
the system is described by some Hamiltonian Hq. Fur- 
thermore, the system may be subjected to an external 
time-dependent perturbation, e. g. an external radiation 
field represented by a Hamiltonian H p (t). Invoking the 
Markov approximation, the motion of the reduced system 
is assumed to be representable in terms of a quantum 
master equation with a generator £ t in Lindblad form 
pfl I2H l22l| (choosing units such that h = 1), 



dt 



p = C t (p) = -i[H + H p {t),p] +V(p). 



(1) 



Since we allow the perturbation to be explicitly time- 
dependent the dynamics does not, in general, provide 
a dynamical semigroup. However, since Ct is in Lind- 
blad form for each fixed t, the time-evolution defined by 
Eq. Q yields a two-parameter family of quantum dy- 
namical maps 23] . 

The commutator in Eq. Q represents the unitary part 
of the motion of the reduced system, while T>{p) is a 
super-operator known as the dissipator of the master 
equation, 



J2jf(AipAT-l{ATAf,p} 



(2) 



The Lindblad operators Af provide the coupling to the 
various decay channels of the open system with corre- 
sponding damping rates jt. We shall assume in the 
following that the master equation Q has been derived 
from an underlying microscopic theory through the weak- 
coupling approximation. This implies that the Lindblad 
operators are obtained as eigen-operators of the unper- 
turbed system Hamiltonian, 



and satisfy A^ 
that 



[H ,Af] =±0JiAf, (3) 
(A~y. The relation © also implies 



[H 07 A+Ar] = [H ,ArA+]=Q. 



(4) 



The damping rates 7^ of the decay modes are determined 
through certain reservoir correlation functions satisfying 
the Kubo-Martin-Schwinger condition which leads to the 
relation 



7i 



exp ( 



(5) 



As mentioned in the introduction the dynamics pro- 
vided by the master equation Q can be represented in 



terms of a piecewise deterministic process ip(t) (PDP) in 
the Hilbert space of the open system in the sense that its 
covariance matrix 

P {t) = E[m))(m\] = I D*pDrpbP,t]\i>)(i>\ w 



satisfies the master equation Here, E denotes the 
expectation value of the stochastic process ip(t) which 
is defined in terms of the corresponding time-dependent 
probability density functional P[ip,i\. 

An appropriate stochastic differential equation for the 
state vector dynamics is given by 



dip(t) 



-iG t (iP{t))dt 



\u7mw 
f A+m 



m)dN+(t). (7) 



The first term on the right-hand side provides the deter- 
ministic evolution periods of the process ip(t). It corre- 
sponds to the deterministic non-linear Schrodinger-type 
equation ip — —iGt{4>), where 



G t {i > )=H^ + l -Y J {l7\\Al 



7+IK + ^l| 2 )V> (8) 



is a non-linear operator with a non-Hermitian Hamilto- 
nian H = H{t) given by 

H = H + H p (t) - 1 £ (j^A+A- + 7+ At At) . (9) 



The quantities dNf(t) in Eq. J7J) are the increments of 
integer-valued random processes N±(t) satisfying 

dN?(t)dN<-(t) = SijSudN^t), (10) 
E[dNf(t)] = jfWAfmfdt, (11) 

where fc, I — ±. According to relation (|10|) these incre- 
ments take on the values or 1. Moreover, if dN^t) 
is equal to 1 for a particular pair (i,k), all other incre- 
ments are equal to zero. Thus, in view of the stochastic 
differential equation Q) the state vector then undergoes 
an instantaneous quantum jump of the form 



\\A>?m\\' 



(12) 



The relation l|ll|) states that the expectation values of 
the increments dNf (t) increase linearly with the time in- 
crement dt and are directly proportional to the damping 
rates 7^ and to the square of the norm of A^tp(t). Since 
the latter depends on time through the deterministic mo- 
tion, the processes N^(t) are, in general, inhomogeneous 
in time. They count the number of quantum jumps of 
type (i, fc), i.e. the number of the jumps (|12fl with Lind- 
blad operator A\. 
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III. DEFINITION OF THE ENTROPY 
PRODUCTION RATE 

A. Entropy flux and numbers of quantum jumps 

An appropriate starting point for the definition of the 
entropy production rate a in an open quantum system is 
provided by the entropy balance equation, 



dS 



(13) 



Here, S is the von Neumann entropy of the open system 
(the Boltzmann constant is set equal to 1), 



S = — tr {plnp} 



(14) 



The quantity J represents the entropy flux, that is the 
entropy which flows per unit of time from the open sys- 
tem to the reservoir, in which case we have J > 0, or from 
the reservoir to the open system, such that J < 0. Thus, 
according to Eq. I|13|) the quantity a is the rate at which 
entropy is produced owing to irreversible processes. 

In order to define the entropy production rate a we 
therefore need an expression for the entropy flux J. To 
this end, we recall that the external reservoir causes the 
wave function to perform instantaneous quantum jumps 
given by (|12[) . According to Eq. © the operator Af in- 
duces a transition of the system state in which the energy 
changes by the amount ±u>i. Thus, Af describes upward 
and A^ downward transitions of the system. Any quan- 
tum jump with jump operator Af {Aj) corresponds to 
the absorption (emission) of a quantum uii of energy from 
(into) the reservoir. Because the reservoir is assumed to 
be in a thermal equilibrium state at temperature T, the 
entropy exchange SS associated with a quantum jump 
can be defined by 



UJi 

5S = —. 

T 



(15) 



To justify this expression one first observes that within 
the Markov approximation of the reduced system dynam- 
ics the reservoir correlation functions decay on a time 
scale which is - by assumption - short compared to the 
relaxation time of the open system. It is therefore rea- 
sonable to assume that the energy exchange in a quan- 
tum jump with operator Af appears as a change ^ftUi of 
the heat energy of the reservoir. We further note that 
the reversible character of the energy exchange through 
a quantum jump derives from the fact that for each jump 
with operator A~ (emission) there also exists the reverse 
jump with operator Af (absorption). 

Thus we see that in each quantum jump with operator 
Af an entropy quantum (I15f) is absorbed from or emit- 
ted into the reservoir. The integer- valued process Nf(t) 
represents the number of quantum jumps with operator 
Af within the time interval [0, t]. It follows that in the 
time interval [0, t], on the average, the net amount 



J2^E[Nr(t)-Nf(t)] 



(16) 



of entropy flows into the reservoir. Consequently, the 
mean entropy flux is obtained from the expression 

• / = E^| E KW-^w]. (17) 



This relation finally leads to the following stochastic ex- 
pression of the entropy production rate of the open sys- 
tem, 

'-f+E?|"[«rw-«fw]. Mi 



By use of this expression the entropy production rate o~(t) 
can easily be determined by a Monte Carlo simulation of 
the stochastic process {7J . To this end, one evaluates the 
von Neumann entropy by means of Eq. ijBJ and records 
the numbers Njf(t) of the various quantum jumps carried 
out during the time interval [0,t]. 



B. Positivity of the entropy production 

We note the following important properties of the en- 
tropy production rate a defined by expression 1)180. First, 
it is immediately clear that a is an extensive quantity: 
When combining statistically independent systems which 
are coupled to the same reservoir, the additivity of the 
von Neumann entropy and the additivity of the numbers 
of quantum jumps is obvious. Second, Eq. (|18|l implies 
that a does not depend on the specific representation of 
the state space of the reduced system, since neither the 
von Neumann entropy nor the number of quantum jumps 
of a given type changes under unitary transformations. 
In particular, the expression l|18|) is invariant under a 
transformation to the interaction picture with respect to 
the system Hamiltonian Hq. Third, the entropy produc- 
tion a as defined in this way is always non-negative, 



a > 0, 



(19) 



which corresponds to the second law of thermodynamics. 

To prove the inequality i|19|) we first write the entropy 
production rate as a functional a = a[p] of the reduced 
density matrix p. Invoking the master equation JQl one 
finds the time-derivative of the von Neumann entropy, 



dS_ 
~dt 



= -tr{2?(p) hip}. 



(20) 



To obtain an appropriate expression for the entropy flux 
we use Eqs. (f 1 Tfl and (f 1 1 1> . as well as the representation 
© of the density matrix to get 



J = J2^tv{ 7 rAtAr p -j+ArAtp} 



(21) 



On the other hand, employing definition @ of the dis- 
sipator and the commutation relations © and (JIJ one 
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finds 

tr {H V(p)} 

= £ 7 rtr U+[Ho,Ar} P+ l[A+A-,H }p\ 

i ^ 

+ 7+tr j^r [ff , A +]p + ^,4+ ff ]p j 
= - ^ tr {tT^ + A~P - X + A- A+p) 

i 

= -TJ. (22) 

The commutation relations Q also imply that the 
Gibbs state pth = exp(—Ho/T)/Z with partition func- 
tion Z represents a zero-mode of the dissipator of the 
master equation, that is 

V(pth) - 0. (23) 

Physically, this means that the Gibbs state is a station- 
ary solution of the master equation without perturbation 
(Hp = 0). We express the Hamiltonian Hq in terms of 
the Gibbs state, 

-H /T = \np th + lnZ, (24) 

which, together with Eq. 1)2 2[l. enables us to write the 
entropy flux as follows, 

J = tr{D(p)lnpth}- (25) 

Combining this with Eq. (|20|l we are finally led to the 
following expression for the entropy production rate in 
the state p, 

a[p] =-tr{V(p) (In p- In p th )}. (26) 

The inequality 1|19[1 can now be shown in two alterna- 
tive ways. One way is to apply a theorem by Lindblad 
[24| to the dynamical semigroup At = exp(2?t) whose 
generator is identical to the dissipator T>(p) of the mas- 
ter equation, and to relate the functional (|26[1 to the 
time-derivative of the relative entropy with respect to the 
Gibbs state pth- Another possibility is to use Lieb's the- 
orem .25] to conclude that the map p 1— > — tr {D(p) In p} 
is a convex functional. Since tr {T>(p) lnp t h} is linear in 
p, it follows that the entropy production rate a = a[p] 
represents a convex functional of the density matrix. 
With the help of this property it is easy to demonstrate 
that the entropy production is non- negative Q. In fact, 
applying the definition of convexity to the functional 
p 1 ► — tv{T>(p) hip} one is led to the following inequality 
which holds for all A e [0, 1], 

-tr {D(\p + [1 - A]p th ) ln(Ap + [1 - A]p th )} 

< -Atr {D(p) ln(p)} - [1 - A]tr {V( Pth ) ln(p th )} . 

Employing the linearity of T>(p) and Eq. Il2.'j[| we get 

-Atr {V( P ) \n(Xp + [1 - A]p th )} < -Atr {V(p) \n(p)} . 

We divide by A and perform the limit A — * to arrive at 

-tr {D(p) \n(p th )} < -tr {V(p) \n(p)} , 

which, by virtue of l|26|) . is equivalent to inequality (11911 . 



IV. DISCUSSION AND EXAMPLES 

A. General properties of the entropy production 

It has been shown in the preceding section that the 
pieccwisc deterministic process (JTJl for the stochastic 
state vector of an open quantum system yields a nat- 
ural definition of the entropy production as a measure of 
the degree of the irreversibility of the dynamics. In this 
definition the entropy flux has been expressed in terms 
of the random numbers of the various types of quantum 
jumps corresponding to the decay channels of the open 
system. 

The given proof of the positivity of the entropy pro- 
duction rate demonstrates not only that the entropy pro- 
duction a is non- negative (see Eq. 11911 ) but also that it 
is a convex functional and that cr(p t h) = is an abso- 
lute minimum of the entropy production (see Eq. H^tty l. 
Hence, in the absence of external fields (H p — 0) the en- 
tropy production vanishes when the equilibrium state is 
reached, in accordance with the second law of thermody- 
namics. 

We remark that for vanishing external fields a coin- 
cides with the negative time derivative of the relative 
entropy with respect to the Gibbs state pth- This is no 
longer true when external fields are present. Contrary to 
the definition of entropy production in terms of the rel- 
ative entropy, the above definition does not assume the 
existence of a stationary density matrix. 

In the case of a non-vanishing external field (H p 7^ 0) 
one expects, in general, that the systems stays away from 
equilibrium and that, therefore, entropy is continuously 
produced. An important physical situation arises (see the 
examples below) if there exists a stationary solution p s of 
the master equation in the interaction picture. Since the 
expression for the entropy production is not affected by 
the transformation to the interaction picture one finds 
that the entropy production in the stationary state is 
given by the constant value 

- S = E^1 E KW-^W], (27) 

i 

where E denotes the expectation value in the stationary 
state. In such a stationary non-equilibrium state entropy 
is thus produced at a constant rate a s . What happens 
physically is that the energy supplied to the system by 
the external field is ultimately transferred as heat energy 
to the reservoir leading to a constant entropy production 
rate. 

It is well-known that the stochastic representation of 
the master equation Q through a stochastic process of 
the form (J2J) is not unique: In general, one can choose, 
in an infinite number of different ways, the Lindblad 
jump operators and the Hamiltonian operator without 
changing the form of the Lindblad generator £ t of the 
master equation However, different such choices do 
change, in general, the stochastic dynamics 0. This fact 
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is well understood in the framework of continuous mea- 
surement theory according to which different choices for 
the jump operators correspond to different measurement 
schemes used to monitor the open system [2(| [271 • Our 
construction of expression (|f 8fl for the entropy produc- 
tion rate presupposes, however, that the jump operators 
are eigen-operators of the system Hamiltonian, sat- 
isfying the commutation relation J3j. It is this property 
which enables one to establish the connection between 
the stochastic expression (|18f) and the functional (|2(jll . 
Physically, it is the eigen-operator property which allows 
one to associate a certain quantum jump of type (i, k) 
with the exchange of a quantum uji of energy with the 
thermal reservoir. 

In the case of a non-degenerate frequency spectrum 
{u>i} the requirement provided by the commutation re- 
lations J2J uniquely fixes the jump operators and, thus, 
the stochastic state vector dynamics. If a certain fre- 
quency uji is degenerate one has the freedom to replace 
the corresponding jump operators Af a (the index a la- 
bels different Lindblad operators belonging to the same 
frequency u>i) by new operators 



where u a p is a unitary matrix. However, in view of 
Eq. (|21[1 it is immediately clear that expression l|18|) 
remains invariant under this replacement, showing the 
uniqueness of the stochastic expression of the entropy 
production rate. 

B. Entropy production rate of a driven qubit 

We conclude the paper by two examples which illus- 
trate several of the features discussed above. The first 
example is two-state system [2^, a qubit consisting of 
an excited state |e) and a ground state \g) with corre- 
sponding energies E e and E g . The qubit is coupled to 
a thermal reservoir and to a single-mode driving field on 
resonance with the transition frequency u> — E e — E g . 
The following discussion will be carried out in the inter- 
action picture. The external field is treated in the ro- 
tating wave approximation such that the corresponding 
interaction Hamiltoninan H p becomes time-independent 
in the interaction picture. 

The stochastic dynamics given by Eq. 10 is defined 
here through two jump operators, given by the raising 
and lowering operators 

A + = \e){g\, A- = \g){e\, (29) 

with stochastic processes AT ± (t) for the numbers of the 
corresponding jumps. The non-Hcrmitian Hamiltonian 
© therefore takes the form 

H = -^(A++A~)-i (j-A+A- + 7 +A-A+) , (30) 



where the Rabi frequency Q is assumed to be real. Equa- 
tion (|27(l then yields the following entropy production 
rate in the stationary state, 

-~tt 3 ^2 + w (31) 

where 7 = 7~ + 7 + . Obviously, cr s is proportional to the 
inversion S3 — p s ee — p s in the stationary state. 



0.3 



0.2 
0.1 




10 20 30 

FIG. 1: The entropy production rate a(t) of a driven qubit. 
Diamonds: Results of a Monte Carlo simulation of the 
stochastic differential equation Q using a sample of 10 5 re- 
alizations with parameters 7 _ /fi = 0.1 and uj/T = 1.0. Con- 
tinuous line: The entropy production rate obtained from the 
solution of the master equation Q. Dashed line: The sta- 
tionary rate <r s according to Eq. 1311 1. 

Figure shows the result of a Monte Carlo simula- 
tion of the PDP defined by Eq. Q) and compares it 
with the solution of the density matrix equation (0 . In 
the stochastic simulation the entropy production rate has 
been determined with the help of Eq. (|18fl . where both 
the von Neumann entropy 5* and the expectation values 
of the numbers N (t) of quantum jumps have been esti- 
mated from a sample of 10 5 realizations of the PDP. The 
initial condition has been chosen to be the ground state 
\g) of the qubit. In the absence of the driving field the en- 
tropy production decreases to zero during the approach 
to the thermal equilibrium state. However, as can be 
seen from the figure, in the presence of an external driv- 
ing the entropy production oscillates and converges to 
the stationary value given by Eq. JsHJ. 

C. A-configuration with a dark state 

As our second example we investigate a two-level atom 
with transition frequency u. Both levels are threefold de- 
generate forming manifolds with total angular momen- 
tum J e = 1 (excited level) and J g — 1 (ground level). 
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We introduce the energy eigenstates \g, m g ) and |e, m e ) 
which arc simultaneously eigenstates of the z-component 
of the atomic angular momentum operator with eigenval- 
ues m g and m e , respectively. In addition to the coupling 
to a thermal reservoir the atom is subjected to a resonant 
driving field which is linearly polarized in the y-direction. 
If we take some initial state in the manifold spanned by 
the states \g, m g = ±1) we then find that the dynamics 
is confined to the subspace spanned by the basis states 



|e,m e = 0), \g,m g = +l) 



\g,m g = -l), (32) 



forming a level scheme known as A-configuration [29| . 

The dynamics of the system will again be studied in 
the interaction picture. It can be described by the four 
jump operators, 



a: 



i 



V2 



!.9,-i)M|, 



An 



1 



V2 



|.g,+l)(e,0|, (33) 



and = (A^y , A\ — (A^y . The corresponding ma- 
trix representation of the non-Hermitian Hamiltonian @ 
in the basis reads 



H = 



7 


9. 


— Q 


-n 


2 





n 





2 



(34) 



A numerical simulation of the dynamics is shown in 
Fig- El In the Monte Carlo simulation of the correspond- 
ing stochastic differential equation Q the quantity a(t) 
has been determined with the help of Eq. I|18(l by record- 
ing the stochastic numbers N^ 2 (t) of upward and down- 
ward transitions in a sample of realizations of the process. 

The stationary entropy production rate may be found 
from Eq. J57J) with the result 



T 



7 P23 



T 



2 7 -+4 7 + + S(7 



+ 



(35) 

This equation shows that a s is proportional to the matrix 
element p 23 = (g, +l\p s \g, — 1) of the stationary density 
matrix p s . We observe that, similarly to the case of the 
driven qubit (see Eq. (|31|) h a s approaches a finite value 
as Q — ► oo. For zero Rabi frequency, Q = 0, the sta- 
tionary entropy production rate vanishes, in accordance 
with our general discussion. Moreover, the entropy pro- 
duction vanishes for infinite temperatures for both the 
driven qubit and the A-configuration. 

By contrast to the driven qubit, however, in the case 
of the A-configuration the limit T — > 0, corresponding 
to a reservoir at zero temperature, leads to a vanishing 
entropy production rate. In fact, in this limit we have 
7~ — ► 70, P23 — ► 1/2 and, therefore, 



lim cr s = lim ^e- w/T 7n^ = 0. 



T^O T 



(36) 



This behaviour can easily be understood if we note that 
the matrix element P23 of the density matrix describes 




FIG. 2: The entropy production rate a(t) of a driven A- 
configuration involving a dark state. Diamonds: Monte Carlo 
simulation of the corresponding stochastic differential equa- 
tion Q with 10 s realizations and parameters 7" /Q, — 1 and 
ui/T = 1.25. Continuous line: The entropy production rate 
obtained from the solution of the master equation. Dashed 
line: The stationary entropy production rate a s given in 
Eq. 1551 . 



the coherence of the state 



|^ d ) = _ +1) + 12,-1)). 



(37) 



This state is known as dark state since it does not couple 
to the driving field on account of a quantum interference 
effect. 

According to Eq. l(3l)|) the stationary entropy produc- 
tion rate is thus proportional to the coherence of the dark 
state in the stationary state p s . At zero temperature, any 
initial state is driven into the dark state under the dy- 
namics. The reason of this fact is that represents 
a zero-mode of the non-Hermitian Hamiltonian 134fl for 
T = 0, that is for 7+ = 0. Equation (|36|) therefore states 
that the entropy production rate vanishes in the dark 
state. This follows immediately from the fact that |^d) 
belongs to the zero temperature ground state manifold. 

Finally, we remark that this system provides an exam- 
ple of a system with a degenerate frequency spectrum, 
mentioned at the end of Sec. IIV Al Both jump operators 
given in (|33|l belong to the same transition frequency ui. 
In a continuous measurement interpretation they corre- 
spond to a measurement of the angular momentum of 
the emitted quanta along the z-direction 30] . The mea- 
surement along any other direction amounts to a trans- 
formation of these jump operators which is of the form 
(|28|l and which does not change the expression for the 
entropy production rate. The latter is therefore seen to 
be independent of the measurement scheme. 
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